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Abstract 

Based on recent work by Adams, I construct a lattice fermion operator that 
fully lifts the staggered flavor degeneracy. The resulting operator is of Wilson 
type but smaller by a factor of 4, better conditioned and contains 3 instead 
of 15 doublers. It is further suggested that this operator may be used as a 
candidate kernel operator to an overlap construction. Prospects for practical 
applications and potential problems of the new discretizations are briefly 
discussed. 

When regularizing a massless continuum Dirac operator on the lattice, 
one is faced with the fermion doubling problem [0, 0]- In addition to the 
single, physical fermion flavor, 15 doubler fermions appear at the edges of 
the Brillouin zone. Traditionally there have been two mutually exclusive 
strategies to ameliorate or solve this problem: One can either remove an 
exact fourfold degeneracy of the naive discretization 0, Q, 111] and reduce the 
doubling problem to 3 doubler species, or, alternatively, one can introduce 
a momentum dependent mass term to lift the degeneracy of the 15 doubler 
species with the physical one jsj. 

In this paper I give an explicit construction of a single flavor fermion 
operator that combines the above mentioned two approaches. It is based on 
recent work by Adams [?J where the two flavor case has been discussed. The 
construction involves adding a momentum dependent mass term as in 
but it starts from the staggered operator [H| rather than the naive fermion 
operator. The resulting operator will be Wilson-like (in particular it will 
break chiral symmetry and require an additive mass renormalization) but 
smaller by a factor of 4 and with a better condition number and only 3 doubler 
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fermions. It is further suggested, that this Wilson-like operator might be a 
suitable overlap kernel operator. 

The starting point of our construction is the massless staggered fermion 
operator \^ 

D st = V ,D„ D, = \ {V, - V$ (1) 

with (r)tj,) X y = (-l)^"^"^ and (V^) xy = U li (x)S x+fijy . This operator obeys 
a remnant chiral symmetry 

{At,e} = (2) 

where e xy = (— iy^i* Xli 8 Xi y. In the spin-flavor interpretation of staggered 
fermions jil, 0L_ e is identified with e = (75 ® £5). 

Following [7J , we first introduce a Wilson-like flavor dependent mass term 

M x = e V5 C (3) 

with (r] 5 ) xy = (viVzVsV^xy = (-l) Xl+X3 8x, y and a Laplacean-like term 
C = (CiC 2 C3C4) sym that is the fully symmetrized product of the = 
I (V^ + V^jj. It is important to note that the spin-flavor structure of Mi is 

M 1 ~(l®&) + 0(a) (4) 

Additionally, Mi has two crucial properties: it is hermitian and commutes 
with e. Both of these follow straightforwardly from the definitions. 

Using theseproperties, we can immediately see that the modified stag- 
gered operator 

D A (m ) = D st + r (1 + Mi) + m (5) 

with the Wilson-like parameter r fulfills a 75-hermiticity like condition 
D A (m )e = eD* A (m ). Consequently, its non-real eigenvalues appear in com- 
plex conjugate pairs which ensures positivity of the determinant for a suitable 
choice of m . Due to its spin-flavor structure (@J, the addition of Mi in (jSJ) 
will spread out the spectrum in the real direction, giving modes a mass term 



according to their approximate flavor chirality (cf. fig. 1(a)). It was demon- 
strated in [?J that this operator is a suitable overlap kernel. The resulting 
overlap operator obeys an index theorem with two fermion flavors 0, flo| . 



1 Note that we have also added a mass term r + mo in order to shift the physical part 
of the spectrum to the correct position. 
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Figure 1: Spectrum of various Wilson- type operators with r = 1 and mo = in the free 
field case on a 32 4 lattice. 



The fact that one is left with two fermion flavors originates in the di- 
mension two of the positive and negative flavor chirality subspaces in four 
space-time dimensions. In order to lift this remaining degeneracy, an addi- 
tional operator is needed, which differentiates between flavors of the same 
chirality. In the flavor Clifford algebra, the natural candidates are the matri- 
ces o^y = i£ v £fi- The a^y commute with £ 5 and can therefore simultaneously 
be diagonalized. Furthermore, has one eigenvalue 1 and one —1 in both 
the positive and negative chirality subspace. Therefore, one can choose a 
common diagonal basis where 

& = diag(l, 1,-1,-1) a ^ = diag(l, -1,-1,1) (6) 

In order to fully lift the flavor degeneracy, we therefore would need an addi- 
tional "mass term" with the flavor structure (1 £g> a^ iu ) + O(a). 

A candidate for such an additional term that has the correct flavor struc- 
ture is given by 

Mf = irf^C^ (7) 

with 

{%u) xy = - (V^) xy = (-l)^ +lXi S x>y for fi<u 

_ 1 (8) 

The spin-flavor structure of can be made explicit by noting that 

rinu = e^VuVv M xv = - (e vll ) = {-l) x » Xv 8 x>y for /i < v (9) 
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We have e^ u ~ (7^7^ <8> 6/£/x) and, up to discretization terms, 
(7/x7v <H> !)• Therefore we see that indeed M^ v ~ (1 ® <x M „) + 0(a). 

It is also straightforward to check that M^ v is hermitian and commutes 
with e. We therefore conclude that M^ v is a valid candidate for a flavor 
dependent mass term and may be used for the construction of a Wilson-type 
operator. 

Let us first look at the operator 

D k (m ) = D st + r(2 + M x + Mf) + m (11) 

which is expected to fully lift the staggered flavor degeneracy. The choice 
of n and v 7^ fi is arbitrary and reflects the ambiguity of defining a single 
flavor in the spin-flavor representation. Note that [M^ , e] = together with 
the hermiticity of M% v implies D^{m)e = eD^(m). Consequently, non-real 
eigenvalues of D^(m) also appear in complex conjugate pairs. 



In fig. 1(b) the spectrum of -D k (0) is plotted for the free case. Note that 



the multiplicity of the real eigenmodes is (4,8,4) as compared to (8,8) for 



the free Da(0) displayed in fig. 1(a), indicating that the flavor degeneracy 
has indeed been lifted. 



In order to investigate the symmetries of the single flavor operator (]11|) . 
we proceed to investigate the transformation properties of under the 



the staggered symmetries [11]. Note, that in addition to breaking the flavor 
chiral symmetry U(l) e , M% v also breaks the discrete shift, axis reversal and 
hypercubic rotational symmetries. Explicitly, we have 

Mr->Jwr-{ ~\ e f ° r e p = ^ y (12) 

under both shift and axis reversal symmetries along the direction p and 

Mf R^R^Mf (13) 



2 Note, that in two dimensions M^ 2 is the unique term required for fully lifting the 
twofold staggered flavor degeneracy (up to a trivial sign). Furthermore, one can obtain 
the flavor chiral mass term M\ as the antisymmetrizcd product of the M^" via 

Mi = ~e aP ^MfM^ (10) 
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under a rotation of the p into the a direction, where the rotation matrix 
is explicitly given by 

{x a for p = p 
-x p for p = a (14) 
x p else 

Note that rotational transformations (Tl3|) can generate an term with 

arbitrary a and f3 ^ a out of any given M 2 V . This suggests that instead of 
adding Mi and a single M% v to the staggered operator as in ( II ip . one might 
instead take a more symmetric linear combination of the different M% v . 
We proceed to investigate the following linear combinations 
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31 n 



M s = —( Sl2 (sis 2 Ml 2 + s 3 s 4 M 2 j 

+ Sl3 ( SlS3 M 2 13 + s 4 s 2 M?) ^ 15 ) 

+ Si4(SlS 4 M 2 14 + S2S3M 2 23 )) 

where the s M = ±1 and s pu = ±1 are arbitrary sign prefactors. One can check 
that the effect of both shift translation and axis inversion in p direction is a 
single sign flip s p — > —s p . Similarly, the effect of a single hypercubic rotation 
is to flip the sign of one single s pu (cf. tabled]). As a result, M s is invariant 
under the following discrete symmetries: 

1. Diagonal shift: 

x^a; + l±2±3±4 (16) 

2. Shifted axis reversal: 

x p ^-x p + p (17) 

3. Double rotation: 

x — > 

r(h»)r(p°) x (is) 

with (/i, z/, p, a) any permutation of (1, 2, 3, 4) 

These symmetries form a subrgoup of the original discrete staggered symme- 
tries shift, axis reversal and hypercubic rotation 

Up to discretization effects, the flavor structure £^ of M s is given by a 
linear combination of the a pu , such that in the common diagonal basis one 
can write 

6 = dia g(1 ,l, -!,-!) ^={^™j OT < 19 > 
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Table 1: Explicit table of sign flips in (fT5|) under the various hypercubic rotations i?( p<T ) 



which indicates, that a single M s alone should fully lift the flavor degeneracy. 
% I therefore propose the following symmetrized Wilson-like operator 

D s (m ) = I)* + r (2 + M s ) + m (20) 

Again, one can demonstrate that D s (m)e = e£)|(m) so that the eigenval- 
ues of D s {m) also appear in complex conjugate pairs. The spectrum of the 



free D s (0) is plotted in fig. 1(c) One can see that it is similar to the one 
of ( TTTT) (cf. fig. 1(b)) with the same multiplicity of real eigenmodes (4,8,4). 
Note however that the spectrum does display a higher degeneracy of the 
eigenmodes which reflects the improved symmetry properties of the free (1201) 
as compared to the free (TTTT) . 

In the interacting theory ( !20l) will receive radiative corrections. The struc- 
ture of these corrections will be restricted by the symmetries (I16H18P (in addi- 
tion to the usual gauge and baryon number symmetries) . These symmetries 



form a subgroup of the discrete staggered symmetries [11] . The counterterm 
structure of the single flavor operator (120)) will be the one of the 2-flavor 
operator (jSJ) plus additional local terms that break the hypercubic rotational 
symmetry but do preserve the double rotation symmetry ( ITBl) . 

The dimension 3 counterterms that can appear are therefore the scalar 
XX j the 2-flavor mass term x^iX an d all of the x^MsX- As i n the case of Wil- 
son fermions, XX necessitates an additive mass renormalization. To leading 
order, the effect of the 2-flavor mass term x^iX w iU be a relative shift of 
the physical and doubler branches of the spectrum as one can see from (jl9j) . 
Mixing among the different xM s % on the other hand will lead to a renormal- 
ization of the flavor structure. This is evident from the fact that the different 
^ do not all commute among each other. The flavor assignment to leading 
order will no more be given by the ^ but instead by a linear combination of 



3 Note that here the two flavors of one chirality receive an opposite mass term while 
both flavors of the other chirality remain unaffected up to discretization effects. 
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themQ The radiative corrections will therefore modify the particular linear 
combination of the four degenerate flavors of the underlying staggered oper- 
ator that will end up in the physical branch of the spectrum. But since the 
details of the flavor assignment are arbitrary in any case and do not carry 
any further physical significance, the renormalized flavor assignment is as 
good as the bare one and there is no need to undo this mixing. Apart from 
this renormalized flavor assignment, the radiative corrections will obviously 
move the relative positions of the physical and doubler branches. This how- 
ever can be absorbed by a further additive mass renormalization as long as 
a clear separation between the physical and doubler branches is maintained 
(which is expected to be the case outside the strong coupling regime). 

Note that radiative corrections are generally expected to be suppressed 
by UV-filtering (and Symanzik improving;) the fermion operator in much the 



same way as it is the case for staggered [12|, [13[ and Wilson fermions [14J, [15 

A further detailed investigation of this issue is beyond the scope of this 
paper but will be an essential ingredient in judging the practical usefulness 
of (JH]). 



We now proceed to use D s (m) as the kernel of an overlap operator [16 
We define a massless overlap operator 

D 1 = p(l + esign(e J D s (-p))) (21) 

that is conjectured to describe a single flavor if p is chosen properly, i.e. in 
a range that extends to p G (0, 2) in the free case. Whether (j2T!) constitutes 
a valid overlap operator remains to be seen. Specifically, the counterterm 



structure and the question of locality and a spectral gap [17| in the hermitian 
kernel operator eD s (—p) require further investigations. 

Beyond the remaining conceptual questions it is interesting to speculate 
about the usefulness of either (T20]) or (TSTj) for lattice QCD calculations. The 
obvious advantage of (I2"U|) as compared to a standard Wilson operator are 
its smaller size (by a factor of 4) and the reduced condition number of D^D 
(another factor of 4 for mo = in the free case). This leads to a reduction in 
storage requirement by a factor of 4 and to a naively estimated speedup factor 
of ~ 16 for conjugate gradient inversions or the construction of the overlap 
operator (12 ip . Note that in contrast to the staggered fermion case these 
improvements originate in the huge reduction of UV modes (by eliminating 12 



'Note however, that all linear combinations of the S,^ do still commute with £5. 
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of the 15 doubler fermion species) and do not come at the price of quadrupling 
the IR modes in the physically relevant branch of the spectrum. 

There are also two obvious disadvantages of ( 120]) as compared to a stan- 
dard Wilson operator. The first one is the appearance of two-hop terms in 
the Wilson-like operator (120]) . These terms could lead to a slower numerical 
implementation when compared to strictly one-hop operators. It is however 
interesting to note, that the additional gauge links needed for the two-hop 
terms are exactly the same gauge links one needs for the construction of a 



clover-term [18| that can be used to 0(a) improve a Wilson operator. One 
therefore expects on the one hand that the speed loss due to 2-hop terms 
would be approximately equal to the speed loss by including a clover term 
(which is typically about 30%) and on the other hand that the inclusion of 
a clover term into (120]) would essentially be free. 

A second concern is that (120]) was constructed using the staggered spin- 
flavor basis and will therefore inherit the nonlocal definition of the spin ma- 
trices. In contrast to staggered fermions however, there will be no significant 
flavor mixing due to the Wilson-like lifting of the flavor degeneracy. Doubler 
states will quickly die out in correlation functions and will not have to be 
disentangled as for staggered fermions. I therefore do not expected that e.g. 
ground state hadron masses or simple matrix elements like decay constants 
or Br will be harder to obtain than with Wilson fermions. It might how- 
ever turn out to be difficult to extract some short distance observables with 
nontrivial spin structure in this formulation. 
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